Entropy exchange in a mixture of ultracold atoms 
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We investigate experimentally the entropy transfer between two distinguishable atomic quantum 
gases at ultralow temperatures. Exploiting a species-selective trapping potential, we are able to 
control the entropy of one target gas in presence of a second auxiliary gas. With this method, we 
drive the target gas into the degenerate regime in conditions of controlled temperature by transferring 
entropy to the auxiliary gas. We envision that our method could be useful both to achieve the low 
entropies required to realize new quantum phases and to measure the temperature of atoms in deep 
optical lattices. We verified the thermalization of the two species in a ID lattice. 
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In recent years an intense research of quantum phases 
common to condensed matter systems and atomic quan- 
tum gases has made remarkable progresses (2]. Some of 
these phases can only be reached provided that the tem- 
perature is suitably low. However, in strongly correlated 
quantum systems, even the temperature measurement 
can be a challenging task. If so, to ascertain whether 
a given quantum phase is accessible, it is convenient to 
focus on the critical value of entropy, rather than tem- 
perature. The advantage is especially clear when the 
strongly correlated regime is reached by sufficiently slow, 
entropy-preserving, transformations of the trapping po- 
tential, as it is often the case for atoms in deep optical lat- 
tices [2|- For these reasons, it is important to determine 
and grasp control of the entropy of degenerate quantum 
gases 0,13) Hi- In this work, we demonstrate the reversible 
and controlled transfer of entropy between the two ultra- 
cold, harmonically trapped Bose gases, which is based 
on the use of a species-selective dipole potential (SSDP), 
i.e., an optical potential experienced exclusively by one 
species (Fig. [T]) [1, Q- In particular, we drive the target 
gas across the threshold for Bose-Einstein condensation, 
by a reversible transfer of entropy to the auxiliary gas. 

The main idea can be understood from textbook ther- 
modynamics. Let us consider two distinguishable gases 
filling an isolated box, exchanging neither particles nor 
energy with the outside, and imagine that only one gas 
(target) is compressed, e.g. through a piston permeable 
to the second gas (auxiliary). The temperature will in- 
crease and, in thermal equilibrium, heat, hence entropy, 
will transfer from the target to the auxiliary uncom- 
pressed gas. In the limit of the auxiliary gas containing 
a large number of particles, it stands as a thermal bath. 
In a more formal way, for an ideal gas of N particles, the 
entropy S is proportional to N log(S/7V), where the num- 
ber of accessible single-particle states S increases with 
the energy density of states and with the average energy. 



i.e., the temperature. In an adiabatic compression of one 
single gas, the reduction of the energy density of states 
is compensated by a temperature raising such that E, 
hence S, remains constant. If we add the uncompressed 
auxiliary gas in thermal contact, the temperature raising 
must be lower: E decreases for the target gas (and in- 
creases for the auxiliary component). In our experiment 
the gases are trapped by adjustable harmonic potentials, 
but the underlying physics is the same. 

To make quantitative predictions we start from the 
entropy of an ideal gas at temperature T in a har- 
monic potential with angular frequency w S = 
kBNth['ig4{z) / g3{z) — log(z)], where iVth denotes the 
number of thermal atoms and the polylogarithmic func- 
tions are defined as gn{z) = X]fe>i •^'^/^"- Above the 
BEC critical temperature Tc, iVth equals the total atom 
number N and the fugacity z is implicitily given by the 
relation N = g3{z){kBT/'huj)^. Below T^, z = 1 and only 
the thermal atoms contribute to the entropy, each with 




FIG. 1: (Color online): schematic of our experimental proce- 
dure. Left: the harmonic magnetic potential is common to 
both gases, auxiliary (red, larger) and target (blue, smaller). 
Right: the species-selective dipole beam compresses the tar- 
get sample and drives it into the degenerate regime. Trapping 
potentials for the auxiliary Rb (dashed) and the target K gas 
(solid) are sketched on the background panels together with 
the K density distributions. 
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a quantity equal to 4fcBC(4)/C(3), Cin-) = .gn(l), so that 

Thus, below Tc the entropy is independent of the total 
number of atoms N. 

The above results refer to a non-interacting Bose gas. 
For high atomic densities, the interatomic interactions 
modify low energy spectrum and must be taken into ac- 
count: a theoretical analysis is given in Ref. 0. We 
evaluate numerically the implicit expression contained 
therein and find that, in our experimental circumstances, 
the interactions increase the ideal gas values of entropy 
by at most 20% for condensate samples. An expression 
for the entropy, easy-to-use and more faithfully approx- 
imating the results of than Eq. H]), is obtained by 
replacing (T/Tc)^ in Eq. ([1]) with the thermal fraction 
1 — /c with fc taken from a semi-ideal model 



_ .3^2/5 
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Here t = T/T^, rj = pL/ksT and the zero-temperature 
chemical potential = {hui /2){lbN as/ ahoY^^ depends 
on as, the scattering length, and aho = yhjrrwj. In 
our experimental conditions, "q ranges from 0.30 to 0.33 
for condensed samples. This approximation is especially 
helpful for high condensate fractions (/c > 0.5), where 
the ideal gas formula Eq. ^ seriously underestimates 
the gas entropy. However, in the explored range fc < 0.3, 
Eq. ID) stands as a reasonably good approximation. For 
thermal samples above the BEC threshold, the interac- 
tion energies are negligible with respect to the tempera- 
ture. This is also the case for interspecies interactions, 
since the auxiliary sample, always thermal for the re- 
ported data, has low density. For this reason, we simply 
add the single-species contributions to obtain the total 
entropy of the mixture, S = Sk + Sm^. 

In our experiment, we load a mixture of *^Rb and 
in a millimetric magnetic trap and sympathetically 
cool the mixture 



The magnetic trap provides a 
harmonic confinement with frequencies {lOx ,ujy,Wz) — 
27r X (24, 297, 297) Hz for K and a factor 1.46 smaller for 
Kb. The geometric mean frequencies uj = (w^-Wj^w^)^/^ 
equal to and we have (wRbi'^K) = 27r x (88, 128) Hz. 
To selectively act on the potential experienced by the K 
atoms (target gas) alone, we use a laser beam tuned to an 
intermediate wavelength between the Di and D2 lines of 
Rb (auxiliary gas) , such that the dipole forces on Rb due 
to these two transitions cancel out The beam, lin- 
early polarized, with a waist of 55 /xm, propagates along 
the horizontal y direction (Fig. [1]). We experimentally 
determine the wavelength value to be 789.85(1) nm, by 
minimizing the efficiency of the Raman-Nath diffraction 
caused by a pulsed standing wave on a Rb condensate 
[l3 |. We also measured the residual potential acting on 
Rb to be T^b/^K = 0.08(1). Such a residual potential 
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FIG. 2: (Color online); BEC fraction of the K sample as 
a function of the K harmonic frequency after the compres- 
sion. Data (circles) are compared to the theoretical predic- 
tions based on the numerical (solid line) and analytical ap- 
proximate (dashed line) solutions of Eq. Q. On the right 
axis, the approximated value S/iNks) = 4(((4)/C(3))(l — /c) 
is shown. The inset displays absorption images of the K sam- 
ple before (a) and after (b) the compression. 



weakly deforms only a small central region of the Rb den- 
sity. From the measured temperature increase caused by 
the SSDP on Rb alone, we estimate a residual effective 
increase of WRb lower than 7%, hereafter neglected. 

On the K sample, the additional confinement induced 
by the SSDP beam is instead a harmonic potential whose 
frequencies add in quadrature to those of the magnetic 
trap. The compression occurs mainly along the weak x 
axis of the magnetic trap (lu^ increases up to a factor 5), 
slightly along the z direction {uJz increases less than 8%), 
and is utterly negligible along the propagation y axis of 
the SSDP beam. Overall, the dipole beam raises the ojk, 
hence the K critical temperature, up to a factor 1.7. 

We raise the power of the SSDP beam with an expo- 
nential ramp lasting 200 ms for maximum compression of 
the trap frequency WK,f = 27r x 216 Hz. The adiabaticity 
is preserved since the ramp is longer than the trapping 
periods of "^^K, i.e., (41,3.4,3.4) ms. Thermal equilibrium 
between the two gases is maintained throughout by mak- 
ing all transformations slow with respect to the inter- 
species collision times, a few ms in typical experimen- 
tal conditions: we verified that the temperatures of the 
two species are always equal, within our statistical uncer- 
tainty. Our observables are the number of thermal atoms 
A'th, the number of condensed atoms Nc and the temper- 
ature T of the two species, that are measured by reso- 
nant absorption imaging after all confining potential are 
removed and the atomic clouds have expanded. For each 
species the entropy is either computed from the measured 
temperature and atom number for samples above BEC, 
or obtained from the condensate fraction below BEC. 

In Fig. [2] we show that, upon selective compression, 
the K sample crosses the BEC threshold, thus its entropy 
manifestly drops. Starting with 9.5 x 10^ Rb atoms and 
2 X 10^ K atoms at 0.410(15) /iK, we measure the K con- 
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densate fraction at different compression ratios. Due to 
the large size of the SSDP beam compared to the K cloud, 
we observe that, in the absence of Rb, the compression 
does not increase the K degeneracy, differently from the 
"dimple" configuration, where a tightly locahzed dipole 
trap is used to bring a single species to BEC 

For a comparison of the observed condensate fraction 
after compression with the theoretical predictions, we 
first calculate the final temperature and then obtain the 
condensate fraction from Eq. The final temperature 
Tf of the mixture after the transformation loka ^k,{ 
satisfies the equation: 

5(r„ UKA, Nk) + S{T,, c^Rb, TVRb) = (3) 

where Ti is the initial temperature. In our experimental 
circumstances, the transformation starts with a thermal 
Rb sample and a K sample at the BEC threshold or be- 
low. Therefore, we solve Eq. ([3]) numerically for Tf. We 
also find an analytical approximate solution for Tf 

Tf = T,(a;K,fM,i){M^[(c^Ka/^K,f)'ee^]/e}'/' (4) 

where ^ = 4C(4)(fcBT,/?iwK,i)V^Rb and W{y) is the 
Lambert's function defined as the inverse of y = xe^ . 
In Fig. [2] we plot the condensate fraction, and the related 
entropy per particle, as a function of the K harmonic fre- 
quency at the end of the compression. The data points 
are well in agreement with the numerical results, that 
differ from the simple analytical solution obtained from 
Eq. ^ and (|4]) by only 10% in our range of compressions. 

We now illustrate how the interspecies exchange of en- 
tropy can be used to explore the entropy-temperature 
diagram of a sample of 1.1(2) x 10^ K atoms, shown in 
Fig. [3] Notice that here we focus on the properties of the 
target K gas, but the presence of the auxiliary Rb gas is 
essential for some transformations. At each data point, 
the entropy of the gas is determined in the following way: 
above Tc, we compute the entropy per particle by the for- 
mula S/Nks = 4 - \og[N{nuj/kBT f] obtained by trun- 
cating the polylogarithmic functions at the first order in 
z] below Tc, we calculate the entropy from the measured 
condensate fraction S/Nks = 4(C(4)/C(3))(1 - /c). 

By combining sympathetic cooling and control of the 
harmonic frequencies, we can follow different trajectories 
in the S-T diagram. In particular, we compare the effi- 
ciency of three trajectories indicated as A, B, C in Fig.[3l 
Trajectory A corresponds to plain sympathetic cooling in 
the magnetic trap. Along B, instead, we raise the SSDP 
beam to maximum power when the temperature of the 
samples reaches 0.6 /iK, and we proceed with evapora- 
tion of Rb and sympathetic cooling of K. When the Rb 
is exhausted, at 0.41 we estinguish the SSDP beam 
intensity, thus decompressing the trapping potential, and 
reach the same final point of the previous trajectory. Fi- 
nally along C, we raise the SSDP when the K sample is 
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FIG. 3: (Color online): trajectories in the entropy- 
temperature diagram of the K sample (see text). For the 
experimental data, the entropy is determined with the for- 
mulas in text. The dashed lines show the ideal entropy for 
a sample of I.l x 10^ K atoms at the harmonic frequencies 
of the magnetic trap (luk — 2-k x 128 Hz) and of the SSDP 
maximally compressed trap {ujk = 2n x 216 Hz). The shaded 
boxes highlight the isothermal SSDP compressions. 



close to BEC at 0.34 fiK and cross threshold by selective 
compression as in Fig.O The end point of all trajectories 
is reached when the Rb sample is either exhausted (B) 
or depleted to the point that further sympathetic cooling 
is inefficient, A^Rb/A^K ~ 2 (A, C). The three described 
paths have similar efficiencies, since the end values of 
entropy are approximately equal. It is important to no- 
tice that, for larger A^Rb/Ax, all cooling processes are 
more efficient and lower S/N values are likely attainable. 
By means of the SSDP compression in presence of the 
auxiliary Rb gas, we perform isothermal transformations. 
Vice versa, with a single species, an adiabatic variation 
of the trapping frequencies corresponds to an isoentropic 
transformation that doesn't increase the gas degeneracy. 

The degree of reversibility of the SSDP transformations 
across the K BEC threshold is investigated by perform- 
ing multiple compression and decompression cycles. The 
SSDP beam intensity is repeatedly ramped up and down 
in an exponential fashion (t = 45 ms) with a period of 
0.43 s. As shown in Fig. [5] the BEC threshold is crossed 
for 5 cycles at maximum compression. The K condensate 
fraction decreases over successive compressions because, 
at each cycle, the number of Rb atoms is reduced and 
the temperature slightly increased. Even in absence of 
compression, our Rb sample experiences an heating rate 
of the order of 0.7 /^K/s which we reduce by an order of 
magnitude by means of a microwave shield [l5| . remov- 
ing trapped Rb atoms with energy E/ks > 5.5 /iK. As a 
consequence, the Rb atom number decreases at a rate of 
2.5(5) X 10^ s~^. Starting with approximately 8.5 x 10^ 
Rb atoms and 2 x 10^ K atoms at 0.4 /iK, we obtain 
a finite K condensate fraction over a time span of 2 s. 
After 5 cycles the entropy exchange still occurs but its 
efficiency is undermined by the lower Rb atom number 
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FIG. 4: (Color online): Cycles of compressions and decom- 
pressions, the K condensate fraction changes in a reversible 
manner as we modulate the K harmonic frequency over time 
from 128 to 216 Hz (the solid line guides the eye). The inset 
shows the relative intensity of the SSDP beam. 



and the temperature before the compression is too high 
to cross the BEC threshold. Our data can still be de- 
scribed by the theoretical analysis summarized in Eq. ([2]) 
and (HI), provided we take into account that successive 
cycles occur from different initial conditions. 

The independent manipulation of two atomic species 
opens also new perspectives for the thermometry of cold 
gases in optical lattices. As a first application, we have 
investigated the thermalization of the two species when 
the SSDP forms a ID lattice for K only The SSDP 
standing wave along the x weak axis of the magnetic 
trap is ramped in 200 ms, kept constant for 20 ms and 
then linearly extinguished in 1 ms. The magnetic trap is 
abruptely switched off immediately afterwards and the 
atoms expand freely for 10 ms (K) and 15 ms (Rb). We 
observe that Rb x and y rms sizes do not depend on the 
SSDP lattice depth and either can be used to extract 
the gas temperature. We also find that the K tempera- 
ture obtained from the y size increases with the lattice 
strength, but it is always equal to Rb within our statisti- 
cal uncertainty of 10%. This shows that, in a ID lattice 
configuration, up to a strength of 20 recoil energy, the 
auxiliary gas allows to read out the system temperature 
and to validate proposed methods to extract the temper- 
ature of a gas in an optical lattice [l3| ■ 

In summary, we have experimentally demonstrated a 
method to exchange entropy between two gases, whereby 
we can precisely reduce the entropy of an ultracold 
atomic sample and drive it across the BEC threshold in 
a reversible manner. Earlier experiments demonstrating 
reversible BEC, or simply gain in phase-space density, 
featured either hydrogen films 17|, a single species gas 
in a "dimple" potential [l3| or multiple spin components 
in the same harmonic potential [l^l ■ Our method is simi- 



components distinguished by their extension. The use of 
two atomic species combined with the SSDP, however, 
offers outstanding flexibility and can be easily extended 
to different configurations, e.g. lattices, and to differ- 
ent mixtures. In particular, the species-selective poten- 
tial works best for mixtures that combine atoms with 
widely spaced Di and D2 transitions and largely differ- 
ent resonant wavelengths. Our method can be applied 
also to reduce the entropy of an ideal degenerate Fermi 
gas, S/N = {Tr^k'^T)/{huj-\/6N) for harmonic trapping. 
In the limit of isothermal transformations, the entropy of 
the target Fermi gas decreases as the inverse of harmonic 
frequency compression ratio Sf/Si — LOi/cuf, less than in 
the case of a condensate Sf/Si = {^i/<^f)^- 

A species-selective optical lattice has important appli- 
cations for the thermometry of an atomic Mott insula- 
tor, filling the need of convenient experimental techniques 
[1^ . Preliminarly, we have verified thermalization to oc- 
cur in a ID SSDP lattice. In addition, the auxiliary com- 
ponent could be used as coolant to dissipate the Mott ex- 
citations [2^. Alternatively, dilute atoms localized by a 
species-selective lattice might be used as disordered scat- 
terers for lattice-insensitive matter waves 
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lar in principle to the "dimple" configuration of Ref. [l3| , 
where the single gas can be thought as consisting of two 
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